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Abstract 

A vortex line is shaped by a zigzag of pinning centers and we study here how far the stretched 
vortex line is able to follow this path. The pinning center is described by an insulating sphere 
of coherence length size such that in its surface the de Gennes boundary condition applies. We 
calculate the free energy density of this system in the framework of the Ginzburg-Landau theory 
and study the critical displacement beyond which the vortex line is detached from the pinning 
center. 

PACS numbers: 74.80.-g Spatially inhomogeneous structures, 74. 25. -q General properties; correlations be- 
tween physical properties in normal and superconducting states, 74. 20. De Phenomenological theories (two- 
fluid, Ginzburg-Landau, etc.) 
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I. INTRODUCTION 

The basic idea of pinning theory is that the vortex is not rigid but adjustable to a local 
distribution of defects^*- In this paper we consider a vortex line pinned by a zigzag of 
pinning centers that are insulating spheres of radius R of the order of the coherence length 
£. The de Gennes boundary condition 4 applies at the pinning sphere surface, considered as 
an insulator-superconductor interface. The straight line connecting these pinning centers is 
a path with abrupt right and left turns, whose zigzag fits into a single plane. As pointed 
out by Yu. N. Ovchinnikov 5 long ago, the investigation of various types of inclusions in 
superconducting materials is of particular interest because of the onset of metastable states. 
We have already considered a similar vortex depinning transition^ to the one considered here 
which is being reviewed from the point of view of displaced defects. 
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FIG. 1: Schematic view of the unit cell middle plane. This view captures the displaced pinning 
spheres of radius R displaced by u along the x direction. The unit cell is orthorhombic with base 
side L and height L z . 

An energetic balance makes the vortex line follows this zigzag path as long as trapping 
by the pinning spheres is advantageous as compared to the increase in length caused by 
the zigzag path. There is a critical path that sets a depinning transition beyond which the 
stretched vortex line is no longer able to follow the zigzag path of the pinning centers. 
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FIG. 2: Nine equivalent unit cells are patched together here and the dark line connecting some 
pinning spheres shows the zigzag path followed by the vortex. 

This critical path is associated to a critical displacement u c that we determine in this 
paper through numerical simulations of the Ginzburg-Landau theory^. From the point of 
view of the Ginzburg-Landau theory, pinning may be caused either by spatial fluctuations 
of the critical temperature, T c (x)i, or by the mean free-path that changes the coefficient in 
front of the gradient term, £(x) 2 |(V — ^A)A\ 2 . The interaction between a vortex line and a 
pinning center has been considered by many authors in the context of the Ginzburg-Landau 
theor y^^ 1 ^ 12 . In this paper we only consider the no magnetic shielding limit, such that 
the field penetrates in the superconductor and there is no Meissner-Ochsenfeld effect. This 
situation can be viewed as a large k limit, k = A/£ being the dimensionless Ginzburg-Landau 
parameter and A is the penetration depth. 

The z-axis is the direction of the magnetic induction, and so, of the applied magnetic 
field, and perpendicular to it is the x-axis. In case of no defects, the vortex is a straight 
line oriented along the z-axis. The defects are equally spaced along the z axis and assume 
alternate displacements u and -u along the x axis. Although the zigzag of defects fit into 
a single plane the problem is genuinely three-dimensional because the defects are coherence 
length size spheres. Consider the position of two defects, (L/2 — u, L/2, Lz/4) and (L/2 + 
u, L/2, Lz/4). All others are obtained from these two forming an infinite lattice in the x,y 
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and z directions. The position of the all others defects can be obtained by translation of 
L along the x or y axis and by Lz along the z axis. The periodicity due to the zigzag of 
defects allow to describe this system through a orthorhombic unit cell of height Lz, and 
of square basis Lx = Ly = L, with two pinning centers inside, at the above positions 
(L/2 — u, L/2, Lz/4) and (L/2 + u,L/2,3Lz/A) with respect to a coordinate frame whose 
origin is at the left-down corner of the unit cell. Figures ^ and El give a pictorial view of 
the geometrical arrangement of pinning centers inside the superconductor which is under 
investigation here. 

The numerical results obtained from simulations of the Ginzburg-Landau theory yield 
the order parameter in each point of space. These results are visualized as iso-surfaces of 
the density, that is of the modulus of the order parameter squared, taken at a fixed value, as 
shown in figures El and 0] Notice that figure 0] corresponds to nine equivalent unit cells, as 
shown in figure 01 patched together. Many other possible pinned vortex configurations are 



FIG. 3: The vortex line remains pinned by the two displaced insulating spheres of radius R = 1.8£ 
in an unit cell with L = Lz = 12£. The picture is the isosurface |A| 2 = 0.2507, which means a 
density value approximately one fourth of its maximum value. The pinning spheres are shown in 
dark gray color. 
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FIG. 4: Nine equivalent unit cells are patched here displaying the vortex lines pinned by the 
spheres, as described in figure [2 

possible in the pinning center arrangement displayed in figure 121 Figure 0] gives one among 
many other, but it is the only one described by the unit cell of figure H3 However these other 
possibilities are not considered here, although the present method is able to treat them by 
searching vortex configurations in larger unit cells containing several pairs of pinning centers. 
We restrict our study here to the minimum deviation of the vortex from the straight line, 
which is well described by the present choice of unit cell. The present study involves some 
parameters related to the unit cell, L z and L, and to the pinning center, u and R. The 
unit cell basis is fixed to L = 12£, and so our goal is to numerically obtain the Helmholtz 
free energy F(u, Lz, R). For each point we obtain a density plot such as shown in figure OJ 
This paper is organized as follows. In section m we discuss out theoretical approach. In 
sections 11111 we show the results obtained through numerical simulations. In section HV1 we 
summarize the main results of the work. 
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II. THEORETICAL APPROACH 



We find numerical solutions of the Ginzburg-Landau theory in the unit cell using a mesh 
grid to describe i t 8 ' 9 ' 10 . The distance between two consecutive mesh points, a, is equal to 
0.5£, consistent to the fact that a must be smaller than the coherence length £, which is the 
minimum physical scale of the Ginzburg-Landau theory. Thus the number of mesh points 
describing the unit cell square basis is P 2 , and so P = 25 since L = a(P — 1). Notice 
that the pinning sphere also imposes a limit on the mesh parameter which cannot be larger 
than the pinning center radius, that is a = 0.5£ < R, otherwise there will be one or no 
mesh point describing the pinning region. In this case the order parameter does not vanish 
inside the defect but just undergoes a drop on its value. The energy density functional 
of the Ginzburg-Landau theory^ is expressed in units of the critical field energy density^, 
H 2 /Att and the order parameter density, A| 2 is dimensionless, varying between and 1. The 
concept of unit cell brings a periodicity to the problem, and so the search for the free energy 
minimum must be done under a constraint, the number of vortices inside the unit cell, v<f)o, 
described by integers: z7 = n x x + n y y + n z z. The magnetic induction, which is the average 
of the local field taken over the unit cell volume, B = f v hd 3 r/V, has a relationship to these 
integers, and in reduced units is B{x) = 2nK^ 2 (n x x/L.L z + n y y/L.L z + n z z/L 2 ). In the 
present paper we only consider a single vortex oriented along the z-axis, hence v — z. 
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The function r(x) is a step-like function used to describe the pinning spheres in this 
approach^. Explicitly we have t(x) = tx(x)t2(x) and 

where is equal to inside and 1 outside the ith sphere. The above explicit representation 
of the t function is necessary for computational reasons and for accuracy we take that 
N = 8. In the limit iV — > oo, the function r tends to the well-known Heaviside function, 
t(x) = & (^~jf^ — l) © ( ^ X ~r 2 ^ ~~ l)- The most significant advantage of the present method, 
is that the free energy functional, Eq. ^ contains the appropriate boundaries conditions to 
the problem. This removes the necessity of solving the theory in two independent regions 
and later applying the Neumman boundary conditions. Besides the present method easily 
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applies to internal regions of any shape, not just spherical, and finds its solution for the 
given normal-superconductor interface^. 




FIG. 5: Kinetic energy and its derivative in the left and right part of the figure, respectively. The 
valleys in the derivative are associated to depinning. The critical displacement is associated to the 
minimum in these valleys. The pinning sphere radius is R = 1.8£ and L = Lz = 12£ 
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FIG. 6: The kinetic energy and its derivative for a system with R = 2.4£, L = Lz = 12£ similar to 
figure [5J Thus the critical displacement depends on the radius of the defect. 

Since we are in the no shielding limit, the vector potential A(x) is determined from the 
condition of magnetic flux quantization inside the unit cell. The vector potential does not 
participate in the minimization process of the free energy density that only takes into account 
the real and imaginary parts of the order parameter. The free energy density contains two 
terms. The first term is the condensation energy density, —t(x) |A| 2 + | |A| 4 , which in case 
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of no vortices (|A| = 1) and no pinning sphere (r = 1) has the value -0.5. The presence of 
a pinning sphere raises the energy since inside it the density vanishes (|A| 2 = 0). And the 
second term, the kinetic energy density, t(x)£ 2 (V — ^A)A . Notice that there is kinetic 
energy in case of no vortices but with a defect. At the insulating-superconducting interface 
r changes from 1 to and this causes a bending of the order parameter, which has some 
kinetic energy cost. 



III. RESULTS 



The critical displacement u c is better observed in the derivative of the kinetic energy 
density^. Figures |3] and |H1 shows the curves dF^n/du versus u for two selected radii, namely 
R = 1.8£, and R = 2.4£, respectively, and in both double hump structure with a 

local minimum between them is seen. For growing displacement u the first minimum of the 
kinetic energy density occurs for no displacement and the second minimum corresponds to 
the depinning transition u c . After the second hump there is a third minimum associated 
to the depinning from the second sphere. Thus detachment of the vortex from a pinning 
sphere makes the kinetic energy density reach a minimum. To understand this consider the 
local maximum that precedes the detachment. For small u the vortex core is superposed to 
the two pinning centers. As long as the vortex is pinned by the two spheres, and follows 
the zigzag path, the order parameter has to adjust around one single common interface 
that involves both the vortex core and the pinning centers. The maximum stretch of the 
vortex line is the longest zigzag path, a configuration that demands the maximum amount 
of circulating current around the vortex. Thus the maximum stretch of the vortex line is 
achieved when the kinetic energy reaches its maximum. 
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Beyond the maximum stretch the vortex unpins from one of the spheres, and the kinetic 
energy decreases for further displacement u. Consequently one expects that the kinetic 
energy undergoes a minimum because following the detachment of the first pinning center 
from the vortex core the kinetic energy must rise again. Notice that in case of the nucleation 
of two independent surfaces from a single one the deflection of the order parameter must 
be taken into account in each one of the new nucleated surfaces. This nucleation leads to 
an extra growth of the kinetic energy because there is a gradient of the order parameter in 
the interfaces. The total interface area has grown because there are two surfaces instead of 
just one after the nucleation. Here the two interfaces are the z = L z /4 pinning sphere and 
the other one formed by the z = 3L z /4 pinning center together with the vortex core. The 
depinning transition is shown in figure 0that shows a sequence of increasing u displacements 
that lead to the depinning of the vortex line from the zigzag path. The same sequence is 




FIG. 7: Sequence of views of the vortex line for several growing displacements u/£ equal to 0.0, 
0.4, 0.8, 1.2, 1.6, 2.0, 2.4, 2.8, and 3.2. The unit cell is cubic, L z = L = 12£, and the pinning center 
is R = 1.8£. The isosurfaces are taken at a density approximately equal for all cases, namely, one 
fourth of its maximum value, |A| 2 ps 0.25. 
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FIG. 8: The same sequence of views of figured is shown here through z = 3L z /4 cuts of that unit 
cell (L = Lz = 12£, R = 1.8£). For the displacement u = 3.2£, the vortex line is still pinned by 
upper the sphere. 

shown again through contour lines along the two planes that cross the pinning spheres in 
half, z = L z /4(figureEl) and 3L 2 /4 (figure EI). 

The behavior of the free energy density versus the displacement u is shown in figure 
for several radii, ranging from R = 1.2£ to 2.8£, and for a cubic unit cell L = Lz = 12£. The 
left graphic in the figure IT21 exhibits the dependency of the critical displacement u c with the 
radius of the spheres. Due to our choice of unit cell parameter L, R is limited to the interval 
1.0£ < R < 3.0£, otherwise are smaller than the vortex core or the pinning spheres will not 
fit in the unit cell. 

The critical displacement is dependent on the height of the unit cell L z , as shown in 
figures fm and in the right side of the figure fT2l for the case that R = 2.0£. Figure ITT1 shows 
a set of free energy density curves versus L z associated to different u displacements in steps 
of 0.2£, ranging from to 6.0£ in the crescent sequence indicated by the arrow. The 
displacement curve is the lowest in energy and the 6.0£ is the largest, in agreement with the 
idea that large pinning centers increase the energy since they bring non-superconducting 
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FIG. 9: The same sequence of views of figure E| is shown here through z = L z /A cuts of that unit 
cell (L = Lz = 12£, R = 1.8£). For the displacement u = 3.21;, the vortex line is detached from the 
lower sphere. 

regions to the unit cell. The right side of Figure shows that unit cells with increasing L z 
allow for larger critical displacements. The zigzag of pinning centers is not so demanding of 
the vortex line for large L z , and so, we expect that large displacements are possible in these 
cases, as confirmed by figures and 



IV. CONCLUSIONS 



The study of the interaction among pinning centers and vortices is crucial to the un- 
derstanding of type II superconductors. For this reason artificially made pinning centers 
are useful and many kinds have been fabricate d 13ll4ll5liei17 though not of the kind consid- 
ered here. In this paper the Ginzburg-Landau theory was numerically solved ona three- 
dimensional mesh to study properties of a vortex line near a zigzag of pinning centers. The 
pinning centers are insulating spheres of coherence length size and the system was described 
by an orthorhombic unit cell containing two displaced defects. We find that there is a criti- 
cal displacement above which the vortex line is detached from a pinning center. Below this 
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FIG. 10: The free energy versus the pinning sphere dislocation for different radii, all obtained for a 
cubic unit cell with L z = L = 12£. The symbols +, o , •, x, n, o, A, y and > mean the radius R/£ 
equal to 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.4, 2.6 and 2.8 respectively. The line of critical displacement 
u c is also shown here. The vortex remains pinned to the zigzag of defects only for u < u c . 

transition the vortex line is pinned by both defects and above by just one. The study of 
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FIG. 11: Several F(u, Lz, R = 2.0£) x Lz curves, associated to different u displacements, are shown 
here. The arrow indicate the ascending values of u, ranging 0.0 to 6.0£ with a step of 0.2£. For 
reference the u = 2.8£ curve has its points magnified here. 
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FIG. 12: In the left side, the dependence of u c with the radius of the sphere is shown {L = L z = 
12£). In the right side, the dependence of u c with the height Lz is presented for a defect of radius 
R = 2.0f . 

coherence length size defects may also be found useful in some other situations. Recently it 
has been shown that the nucleation of such defects is a mechanism to lower the energy of 
the superconducting stated. 
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